ACCELERATING UNIVERSE IN MODIFIED THEORIES 



OF GRAVITY 

B. C. Paul* P.S. Debnathtand S. Ghose 
Physics Department, P. O. North Bengal University 
SiUguri, Dist. : Darjeeling, PIN : 734 013, India 



Abstract 

We study cosmologies in modified theories of gravity considering Lagrangian den- 
sity f{R) which is a polynomial function of scalar curvature (i?) in the Einstein-Hilbert 
action in vacuum. The field equation obtained from the modified action corresponding 
to a Robertson- Walker metric is highly non-linear and not simple enough to obtain 
analytic solution. Consequently we adopt a numerical technique to study the evolu- 
tion of the FRW universe. A number of evolutionary phases of the universe including 
the present accelerating phase are found to exist in the higher derivative theories of 
gravity. The cosmological solutions obtained here are new and interesting. We study 
modified theory of gravity as a toy model to explore the past, the present and predict 
the future evolution. It is found that all the models analyzed here can reproduce the 
current accelerating phase of expansion of the universe. The duration of the present 
accelerating phase is found to depend on the coupling constants of the gravitational 
action. The physical importance of the coupling parameters those considered in the 
action are also discussed. 
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I. INTRODUCTION: 

The last couple of years witnessed a modest progress in our understanding of the ob- 
served universe because of the advent of new cosmological precision tests, capable of pro- 
viding physics at very large redshifts. The luminosity curve of type la supernovae [1], the 
large scale structure [2] and the anisotropy of the CMB [3] favour a spatially fiat universe. 
Recent decade is witnessing a paradigm shift in cosmology from speculative to experimen- 
tal science due to a large number of observational inputs. The has been predicted that the 
present universe is passing through a phase of the cosmic acceleration. It is also believed 
that the universe might have emerged from an inflationary phase in the past. A large 
number of cosmological models were proposed in Einstein's gravity with early inflationary 
scenario in the last three decades which works well. However, the recent prediction that 
the present universe is passing through an accelerated phase of expansion is interesting and 
a proper cause is yet to be understood. It is thought that the cause of the present acceler- 
ation of the universe might be due to dark energy in the universe. However, the concept of 
dark energy in the Einstein's gravity with normal matter or fields cannot be implemented. 
Consequently it is a challenging job in the theoretical physics to frame a theory for cos- 
mological evolution which could address the origin of dark energy also. It is known from 
the cosmological observations that the dark energy content of the universe is about 70% 
to that of the total energy budget of the universe. As mentioned earlier the usual fields 
available in the standard model of the particle physics are not enough to account for the 
huge dark energy reservoir in the universe, a modification to the Einstein's field equation 
either in the gravitational or in the matter sector, therefore, is essential to accommodate 
the present cosmological observations. The issue of dark energy has been taken up in a 
gravitational theory in presence of a cosmological constant [4, 5]. However, the vacuum 
energy density in such a theory remains constant in the course of cosmic evolution and it is 
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also true that there are known contributions into vacuum energy which are several orders 
of magnitude greater than the allowed cosmological values. These observations led us to 
look for an alternative model or a new physics [6-8]. It has been proposed recently that 
gravity itself, if properly modified, could account for the recent cosmic acceleration [9, 10]. 
The standard Einstein's gravity may be modified at low curvature by including the terms 
those are important precisely at low curvature. The simplest possibility is to consider a 
term in the Einstein's-Hilbert action (it may originate from M-theory) [11]. Carroll et al. 
[12] also suggested that such a theory may be suitable to derive cosmological models with 
late accelerating phase. Although a theory with -^-term in the Einstein's gravity accounts 
satisfactorily the present acceleration of the universe, it is realized that inclusion of such 
terms in the Einstein-Hilbert action leads to instabilities [13]. Subsequently, it has been 
shown that further addition of a R'^- term [14] or ln{R) term [15] to the Einstein's gravita- 
tional action leads to a consistent modified theory of gravity which may pass satisfactorily 
solar system tests, and free from instability problem. It is known that the modified gravity 
with a positive power of the curvature scalar (namely, R^-term) [15-18] in the Einstein- 
Hilbert action admits early infiation. The modified gravity with negative powers of the 
curvature in the Einstein-Hilbert action is recently becoming popular as it might effectively 
behave as a dark energy candidate. Consequently the theory might satisfactorily describe 
the recent cosmic acceleration [9-11]. So it is reasonable to explore a theory which could 
accommodate an inflationary scenario of the early universe and an accelerating phase of 
expansion at late time followed by a matter dominated phase. As a result modified theory 
of gravity which contains both positive and negative powers of the curvature scalar (R) 
namely, f{R) = R+ aR^ + (3^ where a and (3 represent coupling constants with ar- 
bitrary constants m and n are considered for exploring cosmological models. It is known 
that the term RJ^ dominates and it permits power law inflation if 1 < m < 2, in the large 
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curvature limit. It may be mentioned here that an inflationary scenario driven by vacuum 
trace anomaly which corresponds to m = 2 and P = was first obtained by Starobinsky 
[17] for describing early inflation. Recently, in the low curvature limit, a number of f{R) 
models have been proposed in order to accommodate universe with late acceleration using 
a modified gravity namely, f{R) = R — ^ with n > [11, 12]. In the metric approach, 
it was shown that the model is not suitable because it does not permit a matter era [19]. 
Recently, it is known that modified gravity namely, f{R) = R + aR"^ is also not cosmo- 
logically viable because it does not permit a consistent scenario accommodating a matter 
dominated era at late time. It is found that instead of matter era, one ends up with a 
radiation era ( a{t) ^ \/t ). On the other hand, R^ model permits matter dominated 
universe {a{t) ~ t^^^) but it fails to connect to a late accelerating phase. In Ref. [19], it 
was shown that the models of the type where Lagrangian density, f{R) = R — ^ with 
ra > and f{R) = aR"^ with m / 1 are not viable for a realistic cosmological scenario as 
they do not permit matter epoch although late acceleration can be realized [20] . Recently, 
modified gravity with power law in R, i.e., f{R) -gravity is examined and found that a 
large class of models including R"^ model does not permit matter dominated universe. 
Capozziello et al. [21] criticized the claim made in Ref. [19]. Tsujikawa [22] derived obser- 
vational signature of f{R) dark energy models that satisfy cosmological and local gravity 
constraints fairly well. The modified /(i?)-gravity is found to be consistent with realistic 
cosmology in some cases [23]. However, no definite physical criteria known so far to select 
a particular kind of theory capable of matching the data at all scales. However, modified 
gravity namely, f{R) ^ e^; or logR may be useful to build a viable cosmological model 
as they permit a matter dominated phase before an accelerating phase of expansion. The 
motivation of the paper is to obtain cosmological solutions considering non-linear terms in 
R in the Einstein-Hibert action. We explore different phases of expansion of the universe 
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from early era to the present and would like to understand the future evolution in the 
framework of higher derivative gravity. The corresponding field equation obtained from the 
above gravitational action is a fourth order differential equation of the scale factor of the 
universe ( a{t) ). As the field equation is highly non-linear and not simple enough to obtain 
an analytic solution we adopt here a numerical technique to solve it. Our approach here is 
similar to that adopted earlier in Ref. [24] which was recently employed by two of us in 
Ref. [25] . In this approach the field equations are first expressed in terms of two functions 
namely, deceleration parameter (g) and Hubble parameter (H) and its derivatives respec- 
tively, which are then solved numerically. 

The plan of the paper is as follows: in sec.II, we set up the relevant field equations in 
the modified theory of gravity, in sec. Ill, cosmological evolutions are predicted in dif- 
ferent models depending upon the coupling parameters of the action adopting numerical 
technique. Finally in sec.IV, we summarize the results obtained. 
II. THE MODEL: 

We consider a gravitational action with non-linear terms in the scalar curvature (R) 
which is given by 



-g d^x (1) 



where SttG = c = 1, ^ is the determinant of the four dimensional metric and R is the scalar 
curvature. Here f{R) is a function of R and its higher power and represents the matter 
Lagrangian. Variation of the action (1) with respect to the metric yields 

where represents the contribution from matter fields scaled by a factor of jtj^ and T^^^ 
denotes the contribution that originates from the curvature to the effective stress energy 
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tensor. Here, T^^ is actually given by 



(3) 



where prime represents the derivative with respect to the Ricci scalar (R) . We are interested 
to study the role of the geometry alone in driving cosmological evolution in this paper, so we 
set Lm = which leads to = in the subsequent sections. The flat Robertson- Walker 

spacetime is given by 



ds' = dt' - a'{t) dr' + r'{d9' + sin' 6 d(j)') 



(4) 



where a{t) is the scale factor of the universe. Using the metric (4) in the field eqs. (2) (see 
also Ref. [24]) we obtain 



a2 - /' 



lif-Rf')-siRf" 



ad' 1 
2- + ^ = -77 
a / 



2^Rf" + Rf"+R^f"'-hf-R/ 



(5) 



(6) 



where an over dot indicates derivative with respect to the cosmic time t. The scalar 
curvature is given by 

(7) 



„ „ [ a a \ 
R = -6 - + ^\ . 



a 

The Ricci scalar R involves second order time derivative of the scale factor a. As the 
eq. (6) contains R terms, one actually has a system of fourth order differential equations 
of scale factor. In the next sections we consider modified gravity with f{R) of the forms 
: (A) f{R) = R + aR^ - , (B) f{R) = R + b ln{R) and (C) f{R) = R + m e^-""^^ to 
explore cosmic evolution. In the above, fj, , a , b , m and n are constants and /x has a 
dimension of R2 [26] i.e. that of (time)~^, a has a dimension of R~^ i.e. (time)^. 



III. COSMOLOGICAL SOLUTIONS : 
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Using eqs. (5) and (6) we obtain 



{HR-R)f -R'f , (8) 



where H = ^ is the Hubble parameter. As both R and H are functions of a and its 
derivatives, eq. (8) is highly non-linear and a differential equation of fourth order in scale 
factor (a(t)). It is not simple to determine analytic solution of the scale factor of the 
universe with cosmic time in closed functional form. Hence adopt numerical technique to 
study the behaviour of the cosmological models based on the parameters of the modified 
gravitational action. For simplicity we consider f{R) of three different forms in the next 
sections. 
Case A : 

In this case we consider modified gravitational action where f{R) is given by 

f{R) = R + aR'-^. (9) 
The corresponding field eq. (8) becomes 

H= r- ^\^-^-£\+a(HR-R) . (10) 



l + 2aR + ^[R^\R R R^j ' 

The above equation is highly non-linear, however, one can obtain asymptotic solutions 
corresponding to different epoch which are 

(i) an exponential expansion a{t) ~ e^°* is permitted when g = — 1 in the early era, 

(ii) an accelerating universe with a{t) ^ &i & later epoch admitting q = The 
variation of q with a{t) in linear evolutionay and exponential phase are shown in figs.l and 
2 respectively. 

Now using deceleration parameter {q) 

aa H 

^ = --2=-^-!' (11) 
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Figure 1: shows the plot of q Vs a{t) in the unit of ln{t). 




0.5 1 1.5 2 



t 

Figure 2: shows the plot of H with t for different D. 
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one can study the evolution of the universe numerically. Since, q is a function of H and its 
derivative, we can re- write the eq. (10) in terms of a second order differential equation in 
q and H to begin with. Since q contains terms with a, one can replace terms with fourth 
order derivative of the scale factor in eq. (10) by q[H]. The functions q and H are, however, 
not independent. The time derivatives in the above equations may now be replaced by the 
derivatives with respect to H using eq. (11). We get the following non-linear differential 
equation 

q" + u{q, H)q'^ + v{q, H)q' + w{q, H) = (12) 



where 



(2g + 4)/ + 216a(g-l)4ij6 

u{q,H) = - 



v{q,H) = - 



(g2-l)[^4_216a(g-l)3/f6]' 

(4g + 7)//^ + 216a{8q + 5)(g - ifH^ 



w{q,H) = 



{q + l)H[n'^ - 216a{q - 1)^H^] ' 
{q - l)[3/x''(2g + 1) + 1296a(g + 1)(q - ifH^ - 36{q - ifH'^] 



{q + l)m - 216a(g -Ifm] 
where the prime indicates a differentiation w.r.t. H and the functions n, v and w depend 
on a and /x^. The above equation, although highly nonlinear, is a second order differential 
equation in q. Here both q and H are time dependent and cannot be solved exactly to 
obtain a known functional form. We solve it numerically in the next section following the 
approach adopted in Ref. [24]. As is a measure of the age of the universe and H is 
a monotonically decreasing function of the cosmic time, eq. (12) may be used to study 
qualitatively the evolution of the universe in terms of q. Since eq. (12) is a second order 
differential equation, to solve it numerically we assume two initial conditions (here, q[H\ 
and q'[H\), for a given value of H. We choose units so that Ho-, the present value of is 
unity and pick up sets of values of q and q' for = 1 (i.e. the present values) from the 
observationally consistent region [5]. We plot q with H for different configuration of the 
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system. As the inverse of H gives us an estimate for the cosmic age, future evolution is 
understood from the region H < 1 and the past from H > 1 in the (g-if)-plane. Since the 
present universe is accelerating, we use a negative q at the present epoch, H = Hq = 1 and 
the universe has entered into this q nagative phase (i,e, acceleration) of expansion, only in 
the recent past. The model may be used to predict also the future course of the evolution 
of the universe. We note the following : 

(i) For a given value of a and fi, the variation of q with H is shown in fig. 3. We choose 
a = 2, fi^ = 12. In the graph time increases from right to left along the horizontal iJ-axis. 
The upper half of the iJ-axis represents decelerating phase and the lower half represents the 
accelerating phase of the universe. It is evidentthat the universe entered into the present 
accelerating phase in the recent past, the rate of acceleration will increase further which 
will attain a maximum thereafter it decreases. There will be slowing down of the cosmic 
acceleration leading to an epoch when the universe expands without acceleration (which is 
transient and occures at H = 0.375) followed by another phase of expansion. In this case 
the expansion of the universe will be accelerating once again. The universe transits from 
decelerating to accelerating phase at H = 1.36. We note that, in this case, the universe 
remains in the accelerating phase once it transits from the decelerating phase. 

(ii) For a given set of values of a and fi we plot q vrs. H for different initial values of 
q'[l] with q[l] = —0.05, when /x^ = 12 and a = 2, the curves are shown in fig. 4. It is 
evident that the universe entered into the accelerating phase in the recent past followed 
by another phase of deceleration. The duration for which the universe transits from the 
present accelerating phase depends on the initial values of g'[l],and the duration increases 
with increasing initial values of ^'[1]. We note the following (i) for q'[l] = 1 the universe 
transits from deceleration to acceleration at H = 1.43 and acceleration to deceleration at 
H = 0.587; (ii) as initial value of q' is increased, the critical point shifts toward left, i.e. 
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Figure 3: shows the plot of q with for /x^ = 12 and a = 2 . Here we choose the initial 
conditions as ^[1]= -0.5, g'[l]=1.655. 



0.4 
0.2 







\^ 

^ \\q'=l 
q'=1.5^ Y 

'\ \ \ 





0.4 0.6 0.8 1 1.2 1.4 

H 



Figure 4: shows the plot of q with H for different initial conditions (/'[I] = 1, ^'[1] = 1.2, 
g'[l]= 1.5 with g[l]= -0.5. 

occurs at a later epoch. 

(iii) We plot q vs. R for different values of a using g[l] = —0.5, g'[l] = 1.2 and /x^ = 12 
which is shown in fig. 5. The figure shows how the nature of the curves depend on the 
initial conditions and that there will be another sign change of g, in the near future. It 
is evident that the universe transits from decelerated phase to accelerated phase only in 
the recent past(i? = 1.4). The universe may transit again to a decelerating phase followed 
by another accelerating phase if a > 0.815 otherwise it will be always in the accelerating 
phase. However, the rate of acceleration depends on the coupling constant a. It is evident 
that as the values of a is increased the corresponding duration of present accelerating 
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Figure 5: shows the plot of q with H for different value of a (a represents a). Here we take 
/Lt^ = 12 and choose the initial conditions as q[l]= -0.5, q''[l]=1.2. 

phase will get reduced, however, the late decelerating phase in future will be enhanced. 
The universe in the remote future once again might enter into an accelerating phase as 
there will be one more sign flip in q. The plot for a = 0.815 is interesting as it decides 
whether the universe will have another phase of acceleration or not. 

(iv) The variation of q with H for different values of coupling constant /x^ is shown in 
fig. 6 for a = 2. It is evident that if inverse Ricci scalar term in the action is absent, then 
it permits a universe which transits from a decelerating phase to an accelerating phase in 
the recent past {H = 1.4) allowing a further sign flip in q leading to a transition from 
accelerating to decelerating phase only. However, for f^'^ ^ 0, an interesting evolutionary 
behaviour of a universe with three sign flips of q leading to a universe from accelerating 
to decelerating followed by a phase decelerating to accelerating and in future deceleration 
to acceleration might happen. However, the rate of acceleration changes with /x. As /x^ 
increases, the period of deceleration to another acceleration in recent future will decrease 
and finally period vanishes at /x^ = 26.46. 

4 

Special Case : We consider here a = i.e., f{R) = R — Equation (12) now takes 
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Figure 6: shows the plot of q with H for different values of /j,'^ (shown by m^) with a = 2 
and choose the initial conditions as g[l]= -0.5, q''[l]=1.2. 

the form 

„// 2g + 4 ,2 (4g + 7) , 3(g-l)(2g + l) , 36(g - l)^^^ _ ^ 

This case was considered earlier by Das et al. [24]. However, the equation (see eq. (12) in 
Ref. [24]) considered by them is not correct as some of the terms are missing. Consequently 
the q vrs. H curve shown in fig. 7 here is found to differ significantly allowing more phases 
of expansion. We note that the universe in the past may have started from a constant 
decelerating phase transits to an accelerating phase thereafter once again because of sign 
flip in q, the universe may transit to a decelerating phase. It has been shown for /^^ = 0.01 
that in future the universe may attain a constant decelerating phase and thereafter there 
will be a sign flip of q once again. As /x^ is increased the corresponding duration of the 
present accelerating phase increases. For smaller values in n the accelerating phase becomes 
shorter. However, there are epoch when g = 1 is attained before and after the present phase 
of acceleration. In all the cases universe will end up with an acceleration having q = —0.5. 
As /x^ increases, the time at which the universe transits from deceleration to acceleration 
occur at an earlier time and the acceleration to deceleration transition occure at late time. 
The universe once again transits from deceleration to acceleration and this will occure 
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Figure 7: shows the plot of q versus H for f{R) = i? — ^ for different value of fi (shown 
by m^). Here we choose the initial conditions as g[l]=-0.5, q''[l]=1.2. 

earlier as /x^ is increased. 
CaseB : 

In this case we consider higher derivative theory much discussed in recent times which 
has the form namely, f{R) = R + h ln(R), to look for a physically relevant cosmological 
model. The relevant differential equation obtained from the field equation is given by 



g + 3 ,2 

~2 T 5 

q^^ — 1 



,^2(,-l) 



(9 + 1) 



6 



1 



0. 



(14) 



The above equation is highly nonlinear, consequently we adopt numerical technique to 
study the evolution of the universe. Here we look for the evolution for different values of the 
coupling parameter h in the action. The curves in fig. 8 are plotted for different coupling 
constants h which are interesting. We note that the universe transits from decelerating 
phase to an accelerating phase in the recent past which will enter into decelerating phase 
once again in future. For 6 < the duration for accelerating phase is found to be shorter 
than that for 6 > 0. However, for h positive the duration of the accelerating phase will 
be longer if h is smaller. In the case of negative 6, the universe is found to land up at 
a maximum possible acceleration at the present epoch, thereafter the rate of acceleration 
will decrease and consequently transits to decelerating phase once again. For 6 > 0, the 
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Figure 8: shows the plot of q versus H for f{R) = R + b InR for different value of b with 
q[l]=-0.5, q'[l]=l. 

maximum rate of expansion will be achieved in near future. 
CaseC : 

In this case we consider gravitational action with f{R) = R + m e^~'^^h Consequently, 

using eq. (8), we obtain the folowing differential equation 

l-6n(q+l)H^ ,2 8q + 5 - 2An(q'^ - l)H^ , 
1 "I ; 9 "I 7 TTT7 5 + 



g+1 



2{q-l){3q + 4) 



- 24n{q - if - 



e^^ -nm 



0. 



(15) 



(g + l)iJ2 "^■"^■^ ^' 3ri^m{q+l)H^ 
The equation is highly non-linear, we adopt numerical technique to solve it as was done in 
earlier sections. We note the following : 

(i) The variation of q with H for for different values of n is shown in the fig.9. It is 
evident that the time of transition of the universe from decelerating to accelerating phase 
depends on n.As n is increased, the time at which universe transits from decelerating 
to accelerating phase is nearer to the present epoch. Also, the duration of the present 
accelerating phase decreases as n decreases. As n increases, the duration of deceleration 
to once again deceleration increases. 

(ii) We plot q vrs. H for different values of m taking n = 0.1 which is shown in fig. 
10. The curves show that the universe at the present epoch entered from decelerating to 
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Figure 9: Here we choose the initial conditions as q'[l]=-0.5, q''[l]=l and m = 4. 
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Figure 10: shows the plot of q versus H for /(i?) = R + m Exp[—nR] for different value of 
m with g[l]=-0.5, g'[l]=l and n = 0.1. 

acceleration and it will switch over to decelerating phase once again in future. The smaller 
values of m leads to a shorter duration of the present accelerating phase. 

From {q — H) curve as drawn in fig. 4, we obtain sufficient data set from the numerical 
plot using Mathematica [27] . Those numerical values may be used to find a closed analytic 
mathematical structure for q and H which can be determined using a polynomial function 
given hy q = TigUiW. We now explore for the mathematical function for q{H) correspond- 
ing to fig. 4 with initial value q'[l] = 1.2. The corresponding approximate analytic function 
may be expressed as 



q[H] = 45.06381 - 515.24544i7 + 2478.58763i72 - 6555.42602i73 + 10556.26293ii'^ 
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Figure 11: Dashing curve is due to polynomial expression and Line curve is from fig. 4 
with q'[l] = 1.2. 

-10835.40545iJ^ + 7153.27419iJ^ - 2949.71637i?^ + 693.08477iJ^ - 70.97952iJ^ (16) 

We compare the curves obtained numerically with that corresponding to the curve fit- 
ted to an analytic function given in eq. (16). The two curves are shown in fig.ll which 
is found exactly superimposed. This comparison holds good only when H is reasonably 
close to one. Thus we outline how one can determine the approximate relation of q with 
H, which may be employed for other curves also. 



IV. DISCUSSION : 

In this paper we obtain cosmological solutions in higher derivative theories of gravity. 
We explore different phases of expansion that are permitted in the higher derivative theories 
of gravity without matter as a toy model. Since the field equations obtained from the action 
are highly non-linear it is not simple to obtain analytic solution in a closed form. We adopt 
numerical technique to understand the present and future evolution of the universe here. 
The relevant field equations corresponding to modified gravitational theories are rewritten 
in terms of two parameters (i) Hubble parameter (H) and (ii) deceleration parameter (q) 
and its derivatives. We plot q vs. H for different values of the coupling constants of the 
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higher derivative gravitational actions namely, (A) a, fx; (B) b and (C) m and n in section 
III to investigate present and future evolution of the universe. It is found that all the 
modified theories of gravity lead to a scenario which accommodate the present accelerating 
phase where the universe transits from a decelerating phase. In some cosmological models 
we note that the universe might once again transit from the present accelerating phase 
to another decelerating phase. In these models one can determine the duration of the 
present accelerating phase in terms of coupling parameters in the gravitational action. In 
figs. (1-10), it is evident that the present accelerating phase of the universe once again 
will change into a decelerating phase. In sections fig. (5-7), new cases are observed where 
yhe universe from the present accelerating phase transits to a decelerating phase and then 
the decelerating phase again changes to an accelerating phase which finally ends up at 
q = —1/2 independent of the initial values of the parameter /it^. Another interesting case is 
evident in figs. (1, 3, 4) where the present rate of acceleration of the universe will decrease 
and attains zero value, where the universe will grow linearly with time as a{t) ~ t for a 
sufficient time duration. Thereafter the universe may once again enter into an accelerating 
phase with an increasing rate of expansion. A detailed study with matter in the modified 
gravity will be taken up elsewhere. 
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